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Stabilizing Quantum States by Constructive Design of 
Open Quantum Dynamics 

Francesco Ticozzi, Sophie G. Schirmer and Xiaoting Wang 



Abstract — Based on recent work on the asymptotic behavior of con- 
trolled quantum Markovian dynamics, we show that any generic quantum 
state can be stabilized by devising constructively a simple Lindblad-GKS 
generator that can achieve global asymptotic stability at the desired 
state. The applicability of such result is demonstrated by designing a 
direct feedback strategy that achieves global stabilization of a qubit state 
encoded in a noise-protected subspace. 



I. Introduction 

The theory of open quantum dynamics has attracted significant 
interest recently due to the fast development of new experimental 
skills to study, and even design, the interaction between a quantum 
system and its environment. In many applications, the dynamics of 
an open system interacting with its reservoir can be described by 
a quantum Markov process. Specifically, let us consider a finite- 
dimensional quantum system described in a Hilbert space H ~ . 
The state of the system is represented by a density operator p on 
H with p > and trace(p) = 1. Density operators form a convex 
set T)(H) C f)(H) with one-dimensional projectors corresponding 
to extreme points (pure states). We denote by 03(H) the set of linear 
operators on H, with $){H) denoting the real subspace of Hermitian 
operators. Throughout the paper we will use \ to denote the adjoint, 
* for the complex conjugate, and [X, Y] = XY - YX, {X, Y} = 
XY + YX to represent the commutator and the anticommutator, 
respectively. The Markovian quantum dynamics is described by the 
following Lindblad-Gorini-Kossakowskii-Sudarshan master equation 
(LME) ffl, |3: 

p = -i[H, p] + V(L, p) = -i[H, p] + Lpl) - i{L f L, p}, (1) 

where H £ $)(H) is the Hamiltonian and L £ 58(H) describes the 
dissipation due to the environment, accounting for the non-unitary 
part of the evolution. 

The mathematical problem we consider is the following: Find a 
pair (H, L) that makes a given density matrix p globally asymp- 
totically stable (GAS) assuming dissipative dynamics of Lindblad 
type (for the standard definition of GAS see 13). This type of 
dynamic stabilization of quantum states is important in quantum 
information processing applications (4), (3- The main result of this 
paper is a constructive procedure to design a pair (H, L) that renders 
an arbitrary mixed state p GAS. This problem of GAS has been 
considered before for pure states, i.e., rank-1 projectors 1101 . t6l . 
(3- Stabilizing pure states is an important task but it is not always 
possible to stabilize a desired pure state. However, there always exists 
a mixed state arbitrarily close to the target state that can be stabilized, 
hence attaining practical stability, that is, stability of a neighborhood 
of the target state. Mixed-state stabilization was considered in (8) 
but without a constructive procedure to find the required dynamics 
(H,L). The L and H we explicitly construct in this Note are in a 
simple form in an eigenbasis of p, tri- and quintdiagonal, respectively, 
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with off-diagonal elements determined by the spectrum of p. The 
diagonal elements a n and h n of L and H are free variables. Further 
analysis shows that p is GAS for most all choices of the diagaonal 
elements a n and h n , and we give explicit conditions on h n to 
guarantee GAS of p. 

In Section [HI] we illustrate by an example how the Lindblad 
generator L obtained this way can be implemented by reservoir 
engineering via direct (Markovian) feedback |9), QO), fH], fl2l . 
where the controlled dynamics has the form of a Wiseman-Milburn 
Markovian Feedback Master Equation (FME)|9|, 1131 : 

p t = -^rJ+ifc+ifFM+M'F), p t ]+V(M-iF, p t )+V{L Q , p t ). 

(2) 

The drift Hamiltonian H £ $)(H), the measurement operator M £ 
23(H) and the noise operator Lq are assumed to be given, while the 
open-loop and the feedback Hamiltonians Hc,F £ f)(H) are our 
control parameters. 

II. Main Results 
A. Preliminaries and design assumptions 

We recall two results on subspace stabilization that have been 
derived in (3, (3, which will be used in the rest of the paper. 

Theorem 1 (Feedback-attractive subspaces fiT}/): Let Hi = Hs © 
Hr with ns being the orthogonal projection on Hs- If we can freely 
choose both the open-loop and the feedback Hamiltonian then, for any 
measurement operator M, there exist a feedback Hamiltonian F and 
a Hamiltonian compensation H c that make the subsystem supported 
by Hs invariant and attractive for the FME ((3 iff 



^0. 



(3) 



A constructive proof is provided in |7 |. We will also make use of the 
following characterization of the global attractivity of a state: 

Theorem 2 (Uniqueness equivalent to GAS [8)1): A steady state of 
is GAS if and only if it is unique. 

Finally, we require a few basic observations and a simple lemma. 
Consider an orthogonal decomposition H = Hs © Hr with 
d s = dim(Hs), d r = dim(Hfl), and let {|0f )}f =1 , )}^ 
be orthonormal bases for Hs, Hr, respectively. The basis 

{|Vm>} = {|^>}fc 1 U{|0f>}fl 1 , 

induces a block decomposition for matrices representing operators 
acting on H: 

X= ^— . (4) 



r x s 


x P ' 


[ x Q 


Xr \ 



and we have the following: 

Lemma 1: Assume p = [ P q ° ] , divided in blocks accordingly 

to some orthogonal Hilbert space decomposition H = Hs © Hr. 
Then p is invariant for the dynamics ([TJ if and only if: 



= 



\{L\Ls, ps} + LpppLp 



(5a) 



p + L' q Lr) 



= 



i[H s ,ps] + L s psL s 

- H L h L Q'Ps} 

— i(HppR — psFlp) + LspsLg — \ps(LgLi 
+ L pPr L r -\{l\L p + L%Lr)pr (5b) 

i\H R , p R ] + L R p R L R - \{L R L R , p R } + Lqp s V q 
-\{L p Lp, P r} (5c) 

Proof: By direct computation of the generator l[T) one finds its 
S, P and i?-blocks to be the l.h.s. of (f5at-(f5cl> . respectively. A given 
state p is stationary if and only if C(p) = 0, and hence if and only 
if its blocks are all zero. ■ 
We know from |6|, |7| that for p to be invariant, its support H p must 
be an invariant subspace, and using the constructive procedure used in 



2 



the proof of Theorem Q] we can construct an block-upper-triangular 
L that stabilizes the H p subspace. There are many possible choices 
to do that, e.g. 



L = 



' L p 


L p ,p 





Lp.R 



with blocks 







-tp,R 
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with ii, £%,... 7^ 0. Therefore, we can focus on the dynamics 
restricted to the invariant support H p , and restrict our attention here 
to full-rank states p = diag(pi, . . . ,Pn) with pi,... ,pn > 0. To 
develop a constructive procedure to build a stabilizing pair (H, L) 
with a simple structure we make a series of assumptions and design 
choices: 

Assumption 1. The spectrum of p is non-degenerate (generic case). 
A state with non-degerate spectrum can be chosen arbitrarily close 
to any state. Without loss of generality, we can choose a basis such 
that p is diagonal p = diag(pi, . . . ,pjv) with pi > . . . > pjv > 0. 
This assumption is instrumental for the construction of L but can 
actually be relaxed, as we remark after the Theorem [3] Consider the 
decomposition H = Hs ffi Hr, with dim('Hp) = 1, such that the 
corresponding block decomposition for p is p 
accordingly H, L. 

Assumption 2. ps satisfies 



P o P ° R J > and divide 



i[H s ,p s ] + L s psLl - \{L\L s , ps} = 0. 



(6) 



Condition $6$ is clearly not satisfied in general, but it will be ensured 
at each step of the iterative procedure outlined in the next Section. 
Given this, conditions f5aV(|5ct can be rewritten as: 



= LqpsIs q - p R L P L P 

= L P p R L P - ±{LlLQ,ps}. 



(7) 



Since we assumed Hr to be one-dimensional, Lp and Lq are both 
n — 1 dimensional vectors. Call rip = LpL P /(L P Lp), Hq = 
LqLq/(LqLq). Then the second equation in (7) reads: 



pp(L P Lp)n P = (LqLU)^(Uqps + p s n Q ) 



(8) 



Assumption 3. Choose Hq to be the orthogonal projector on the 
eigenspace corresponding to the smallest eigenvalue of p$. 
Choose an ordered spectral decomposition with rank-one projectors, 
p = Y^fJi Pi^s,% + PrRr, where > p i+1 . As p is in diagonal 
form this means that we are choosing: 







We thus get: 



Pr(V p L p )Up 



.. l Q ). 
( l qLq)pn-iU-q, 



(9) 



which can be satisfied if and only if lip = TIq and Pr(L p Lp) = 
(LqLq)pn-i- We can choose L P Lp — pjv-i, obtaining 
(LqLq) = pn- Hence £q can be in particular chosen to be real, 
£q = <JPn- We have thus constructed a Lindblad term of the form: 














Ls 






^PN-1 


o ■ 


■■ 




Lr . 



As pR is a scalar, we can rewrite J5bl > as: 

-i{pRls +ps)H P + K = 0, 



(10) 



where K = -±(L P lS R + L\L p + L t Q L R )p R - \p s (LlL P + 
LqLr) + LspsL^. Once p s , Pr, L s , Lr, L p , Lq are fixed, dTOt 
is a linear system in Hp which always admits a unique solution 
Hp = i(ppls + ps)~ 1 K, which is a necessary condition for the 
invariance of p. We are now in a position to present an inductive 
procedure for constructing stabilizing generators. 



B. Constructive algorithm and proof of uniqueness 

We start by constructing a stabilizing generator for the two-level 
case. Let p = [ P q P2 ] with p\ > P2 > 0. Note that in the reasoning 
of the previous Section it is not necessary to impose tr(p) = 1. 
We will make use of this fact in extending the procedure to the n- 
dimensional case. 

Given the previous observations, we can render the given p 
invariant for the dynamics if we can enact dissipation driven by a 
Lindblad operator and a Hamiltonian of the form 



Oi 

JP2 



02 



(ID 



and a H such that its off-diagonal elements satisfy d 1 0b . 

In the N -level case let p = diag(pi, . . . ,pjv). We can iterate our 
procedure by induction on the dimension of Hs- We have just found 
the 2x2 upper-left blocks of L and H such that p^ = diag(pi,p2) 
is stable and attractive for the dynamics driven by the reduced 
matrices. Assume that we have some m x m upper-left blocks 
of L and H such that p^ 71 ' = diag(pi, . . . , p m ) is invariant for 
the reduced dynamics. This is exactly Assumption 2 above. Let 



H 



(m) 



s = s upp(E}liP3 n j) 
want to stabilize p < - m+1 '> = di 
( m ) m-iy( m ) 



u 



(m) 



= supp(p m+ in m+ i). If we 
. . . ,p m +\) for the dynamics 

- (m) T (m) 



restricted to T-Cg' 1 ' ©"H^"', we can then proceed building Lq , L p 
using Design Choice 1 above. Design Choice 2 lets us compute the 
off-diagonal terms of the Hamiltonian, while we can pick Hr to 
assume any value, since it does not enter the procedure. 

By iterating until m = N, we obtain a tridiagonal matrix L with 
i«,n+l = y/P n , L n +l,n = VP n + 1 and L n ,n = O-n, i.e., 



L = 







y/Pi 

Vp 3 





a-3 



VP N -i 




y/P. 



OJV-1 

VPn 

and H becomes a quintdiagonal Hermitian matrix, i.e. 

H nn ■ — h n j 



N-l 

a N 



(12) 



H n ,n+± — a. 



2 Vp n + Vp. 



n+1 



tln,n+2 — rl n+2 ,n ~ ~7:Pn+l—p j — , 

2 VP n + VP n +2 

(13) 

Theorem 3: If (H, L) are chosen as in d 1 244 1 3 1 then p — 
diag(pi, . . . ,Pn) is a stationary state of the LME p(t) — —i[H,p] + 
T>(L, p). p is GAS for most choices of the diagonal elements a and 
h of L and H, respectively; in particular there exists Mq > so that 
p is GAS for h = (M, 0, . . . , 0) for all M > M . 
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Proof: Given the generator, we can verify by direct calculation 
that p is a steady state. Setting B = T>(L, p), direct calculation shows 
that B,„.n — except for 

1 2 
B n ,n+1 = B n +l,n = -^(V^n _ VPn + l) ( d «\^„ + + 

1 2 
B n ,„+2 = -B n + 2 ,„ = -^Pn + ^V^n _ VP n + l) I 

and setting A = i[iJ", p] shows that the Hamiltonian term exactly 
cancels the non-zero elements of B, i.e., —A + B — 0. Thus p is a 
steady state of the system. 

To show how to make p the unique, and hence attractive, stationary 
state, assume that p' is another stationary state in the support of p. Let 
X = {X\X = xp + yp', x,y G R}. Then any state in X is 
also stationary. Since X is unbounded while D(H) is compact, there 
must be a stationary state pi at the boundary of T>(H), i.e., with 
rank strictly less than p. Then the support of pi, Hi — supp(pi), 
must be invariant |7 |, and H, L must exhibit the following block- 
decompositions with respect to the orthogonal decomposition H = 
Hi 9% : 



H 



Hn Hi 



H 



H22 



Ln 




L12 
L22 



with H12 — ~i\h\ x Li2. If L12 / then it is straightforward to 
show that Tr(Il2p), where II2 the orthogonal projection onto H2, is 
strictly decreasing [6|, and thus H2 is not invariant. Since p is also 
stationary, and is a full-rank state, this is not possible and we must 
therefore have L12 = and thus H12 = 0. This implies that H and 
L are block-diagonal. Both Hi and H2 must contain at least one 
eigenvector of L, say vi and V2- Orthogonality of the subspaces Hi 
and H2 along with the block-decomposition of H and L imply that 
any pair of vectors Vi G Hi, V2 G H2 must satisfy: 



vlLvi = vlHvi = vlvi = 0. 



(14) 



These are necessary condition the existence of another stationary 
state, and hence for making p not attractive. 

From the results in the appendix (Theorem[7]l, a tridiagonal L in the 
form l !12t has n distinct eigenvalues corresponding to n eigenvectors 
Vk with real entries and the first all different from zero. Without 
loss of generality, we can assume that the first element of each 
(unnormalized) eigenvector Vj equals 1 for all j. If the eigenvectors 
are mutually non-orthogonal, i.e., v^ve 7^ for all k, I then the third 
equality in d 14b is automatically violated, and hence p must be the 
unique stationary state. This condition will almost always be satisfied 
in practice, and it is easy to check that it always true when N — 2 
(see appendix). However, even if L has orthogonal eigenvectors we 
can use our freedom of choice in the diagonal elements h n n of H 
to render p the unique stationary state. 

Assume ujffc = for some pair of eigenvectors of L. Let 
Ho be the Hamiltonian corresponding to h = and Mo = 
maxj,fc \ VjHoVk\- Choose M > Mo and let H be the Hamiltonian 
corresponding to h = (M, 0, . . .), H = Ho + diag(M, 0, . . . , 0). 
Recalling that the first component of Vj is 1 for all j, shows that 

v]Hv k = M + v]H v k > 0, 

and therefore the second equality in d 1 4b is violated, and p is the 
unique stationary state. ■ 
Remark: Theorem [3}j7] further shows that Assumption 1 on the 
spectrum of p can be relaxed, since the fact that the spectrum is non- 
degenerate plays no role in the proof. The construction is effective 
for any full rank state on the desired support. 



III. Feedback stabilization of encoded qubit states 

Consider a system whose evolution is governed by the FME l(2j 
with Ho — —H' c , M = M' , and L admitting an eigenspace Hs 
of dimension 2 for some eigenvalue As- Assume we can switch off 
the measurement and the feedback Hamiltonian, F — M — 0. With 
this choice, ((2} admits a two-dimensional noiseless (or decoherence- 
free) subspace |14|, |15|, |6|, which can be effectively used to 
encode a quantum bit protected from noise. We now face the problem 
of initializing the quantum state inside the DFS: we thus wish to 
construct He, F and M, such that a given state p of the encoded 
qubit is GAS on the full Hilbert space. Setting H = Hs © Wr 
and choosing an appropriate basis for Hs, the encoded state to be 
stabilized takes the form 



P = 



The dynamical generators can be partitioned accordingly, 



PS 


0" 


, PS = 


Pi 


0" 











P2_ 



L = 



F 





F S 
Ft 



L P 
Lr 

F P ~ 

Fr 



H c = 



H s 
Hi 



M 



Ms 
Mi 



Hp 
H R 
M P 

Mr 



where I2 is the 2x2 identity matrix. We compensate the feedback- 
correction to the Hamiltonian by choosing He = H' c — \ (FM + 
M ] F). For L P ^ we use the constructive algorithm described 
above to render Hs attractive by choosing Fp = —iMp, and 
constructing a H' B for H' c so that no invariant state has support in 
Hr (see Theorem 12 in |7|) and by imposing H' P = 0. For Lp — 
we need to choose an observable such that Mp 7^ 0. We are thus 
left with freedom on the choice of Hs , Fs , which can be now used 
to stabilize the desired ps in the controlled invariant subspace Hs- 
Denote the elements of the upper-left 2x2 blocks as 



M s = 



'Mi 


M 3 ' 


, F S = 


~Fi 


F, 


Ml 


M 2 _ 




F 2 



(15) 



Set F s = -^fHv^ 2 - v^i)- If M 3 + °> settin g fe M := 
2M 3 /(yp 1 + ^p 2 ) shows that M 3 = ^(VPi + VP 2 )' Le - U P to 
the multiplicative constant Ism the two block matrices are such that 
Ms — iFs is of the form d!2t . Theorem [5] applies and the desired 
state is GAS. Notice that if M3 = 0, and M has non-degenerate 
spectrum, it is easy to find another state p' , arbitrarily close to p, 
such that, in the basis in which p' is diagonal, we have A/3 7^ 0. 
Thus we can attain practical stabilization of any state of the encoded 
qubit. 

IV. Conclusions and outlook 

Efficient quantum state preparation is crucial to most of the phys- 
ical implementations of quantum information technologies. Here we 
have shown how quantum noise can be designed to stabilize arbitrary 
quantum states. The main interest in such a result is motivated by 
direct feedback design for applications in quantum optical and opto- 
mechanical systems and quantum information processing applica- 
tions. As an example we have demonstrated how to devise the (open- 
and closed- loop) control Hamiltonians in order to asymptotically 
stabilize a state of a qubit encoded in a noiseless subspace of a 
larger system. Further study is under way to address similar problems 
for multiple qubits, in the presence of structural constraints on the 
measurements, control and feedback operators, and to optimize the 
speed of convergence to the target state. 
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Appendix 

We collect here results about orthogonal polynomials and tridiago- 
nal matrices that are instrumental to the proof of Theorem [3] Details 
and (missing) proofs can be found in |16|, |17], |18|. 

Definition 1: An orthogonal polynomial sequence {P n (x)}^ =1 
over [a,b] C R is an infinite sequence of real polynomials such 
that (Pj,Pk) = for any j ^ k under some L 2 inner product 
(Pj,Pk) := J Pj(x)Pk{x)w(x)dx with a weighing function w(x). 

Theorem 4 (Favard U6V ): A monic polynomial sequence 
{P n (x)} is an orthogonal polynomial sequence if it satisfies a 
three-term recurrence relation 



Pn+i(x) = (x-b n )P n (x) 



i(x),n>0 (16) 



with {b n }, {c n } sequences of real numbers and c n > 0. 

Theorem 5 ( M7V ): If {P n (x)} is a sequence of orthogonal poly- 
nomials then each P n (x) has n distinct roots. 

Theorem 6: A real-symmetric tridiagonal TV x TV matrix with non- 
zero entries on the first sub/superdiagonal has TV distinct eigenvalues. 

Proof: It suffices to note that the monic polynomial sequence 
{/„} with / (A) = 1 and 



/„(A) = det 



A 



cei 



Pi 











Pi 




A — a 2 P2 



Pn-l 
A — a„ 



••• p„-i 

for n — 1, 2, ... TV satisfies the recurrence relation 

/ n+1 (A) = (A - a n+1 )/»(A) -/£/»_i(A) 



(17) 



(18) 



and therefore is an orthogonal polynomial sequence by Favard's 
theorem provided P„ > 0. Thus f„(X) has n distinct roots for all n, 
and the matrix has TV distinct eigenvalues. ■ 
The expression for the corresponding (unnormalized) eigenvectors 
Vk can then be obtained by straightforward calculation: 

/j-i(Afc) 



(19) 



Pi--- Pi-i 

Theorem 7: Let T be a real tridiagonal matrix. If T n , n +i = P n > 
and T n +i tn = 7„ > then T has TV distinct eigenvalues A n . 

Proof: There always exists a diagonal matrix D = diag(d„) 
such that S — D~ 1 TD is a symmetric tridiagonal matrix: the 
diagonal elements must satisfy d n+1 = d n j n //3 n for all n, which 
shows that D is uniquely determined up to global factor. We fix 
d n = 1. The off-diagonal elements of S are Sn,n-\-l — Sn+l,n — 
V Pnjn- The matrix S is real-symmetrix and hence S = VEV, 
with E = diag(A„). By the previous theorem the eigenvalues A n 
are real and distinct provided /? n 7n > for all n. V is a real- 
orthogonal matrix whose columns are the normalized eigenvectors of 
S, Vk = C^ l Vk with Vjk as defined in d 1 9b and C\ = \\v~k\\- Since 
T = DSD' 1 = D(VEV T )D~ 1 = (DV)E(DV)~ 1 , T has the 
same eigenvalues as S, with eigenvectors Wk = Dvk- ■ 

While it is not strictly needed for the result in this work, it is 
worth noticing that for almost all choices of diagonal entries in 
T, the eigenvectors will be mutually non-orthogonal unless T is 
symmetric. For fixed off-diagonal elements Pu and fk, the eigen- 
values Afc and corresponding eigenvectors Vk of S can be expressed 
explicitly in terms of the diagonal entries a. Since /;v(A, a) is an 
TVth order polynomial in A with TV distinct roots, by the implicit 
function theorem, all roots \k for /jv(A,a) = can be expressed 
locally as continuously-differentiable functions Afc(a) in some open 



neighborhood J\f a of a, and similarly for the eigenvectors Vk- To 
determine whether T(a) has orthogonal eigenvectors consider the 
functions 



F k i{a) = ^2 d£&£-i/n-i(A fe (a))/„_i(Ai(a)) 



(20) 



which correspond to the inner products of the unnormalized eigen- 
vectors of T with b n — (Pi . . . Pnji ■ ■ ■ 7n) _1 - By orthogonality of 
the eigenvectors of S we have 

JV-l 

B U = 6 n-i/n-i(A fc (a))/ n _ 1 (A < (a)) (21) 

n=l 

for all ex. But Fki and Gki are polynomials in A and a and differ 
only by the coefficiently d n . Thus, we will have Fke(oc) ^ Gkiipt) 
for most a unless all d n are equal, which implies p n = 7„ for 
all n and T symmetric. (In our case this would only happen when 
all populations p n are equal.) By continuity of Fke(a) we can 
conclude that in any neighborhood Af ao of a bad point ceo with 
Fke(oto) = 0, there is an open (and thus positive-measure) subset of 
points with Fm(cx) 7^ 0, provided Fkt(cx) is not constant on Af ao . 
This argument can be made stronger and more rigorous if we assume 
9 g"^°^ (cto) 7^ for some a m , i.e., that not all partial derivatives 
at the bad point ao vanish, in which case we can show that Fm(ol) 
vanishes only on a measure-zero subset of M ao using the implicit 
function theorem. 
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